Tailbeat perturbations improve swimming efficiency by reducing the phase lag between body motion and the resulting fluid response

Abstract Understanding how animals swim efficiently and generate high thrust in complex fluid environments is of considerable interest to researchers in various fields, including biology, physics, and engineering. However, the influence of often-overlooked perturbations on swimming fish remains largely unexplored. Here, we investigate the propulsion generated by oscillating tailbeats with superimposed rhythmic perturbations of high frequency and low amplitude. We reveal, using a combination of experiments in a biomimetic fish-like robotic platform, computational fluid dynamics simulations, and theoretical analysis, that rhythmic perturbations can significantly increase both swimming efficiency and thrust production. The introduction of perturbations increases pressure-induced thrust, while reduced phase lag between body motion and the subsequent fluid dynamics response improves swimming efficiency. Moreover, our findings suggest that beneficial perturbations are sensitive to kinematic parameters, resolving previous conflicts regarding the effects of such perturbations. Our results highlight the potential benefits of introducing perturbations in propulsion generators, providing potential hypotheses for living systems and inspiring the design of artificial flapping-based propulsion systems.

Being so fundamental, how fish optimize swimming kinematics to enhance swimming efficiency has been widely explored over the past decades (23).A prevailing understanding from biological (24)(25)(26)(27) and physical (28,29) investigations has revealed that swimming in a Strouhal number (defined as St = 2fA/U, where f is the swimming frequency, A denotes tailbeats amplitude, and U refers to the swimming velocity) of 0.2-0.4allows individuals to achieve optimal swimming efficiency.Besides Strouhal number, other studies have shown that swimming efficiency also benefits from other kinematic properties, such as the amplitude of tailbeat (30), the reduced frequency (normalized swimming frequency) (31), and the profile of oscillating waveforms (e.g.nonsinusoidal (32,33) and quasisymmetric oscillations (34,35)).Saadat et al. (30) revealed that a Strouhal number of 0.2-0.4 in combination with an amplitude of approximately 20% of the body length can lead to efficient propulsion.Floryan et al. (31) found, within a specific Strouhal number, the smaller the reduced frequency, the higher the propulsion efficiency.
Most current studies on swimming kinematics and performance assume that the swimmer's movements can be described using various idealized sinusoidal oscillations.However, in nature, swimmers do not use such "clean" kinematic oscillations to move forward.With fast Fourier transform algorithm, natural kinematic oscillations can be resolved into basic and high-frequency components.In most studies, the high-frequency components, which are typically low in amplitude, are often disregarded in simplified body wave models due to their similarity to noise (36).As a result, few studies have investigated the role of these highfrequency and low-amplitude perturbations in kinematics (37,38).Lehn et al. (37) found that combining high-frequency and lowamplitude perturbations to flexible flapping foils can increase thrust and efficiency.However, it is unclear whether this improvement is mainly due to the flexibility of the foil or the perturbations, as certain stiffness of the foil can also significantly improve efficiency (15)(16)(17).Furthermore, some studies that examined the performance of a flapping rigid foil found that perturbations can enhance thrust but not the efficiency of swimming (38), thus calling into question the mechanism that allows for increased efficiency with perturbations.In addition, previous studies have primarily investigated how perturbation frequency impacts swimming efficiency, while there is a lack of research exploring the effect of perturbation amplitude, to thrust production and efficiency, as well as the combined impact of both frequency and amplitude.Overall, it is still largely unknown whether a swimmer with perturbations can improve thrust and efficiency simultaneously, and if so, what kind of perturbations are beneficial and how they lead to improvement.In this work, we combined experimental, numerical, and theoretical studies to systematically explore the benefits of perturbations over the swimmer's performance (Fig. S1).

Robotic fish experiments
To investigate the general function of small kinematic perturbations to tailbeats during fish-like swimming, we consider perturbations with high frequency and low amplitude compared to a basic sinusoidal wave motion performed by the tail (Fig. 1A and B).
Three modes are defined: basic mode (BM) with a sinusoidal wave representing rhythmic tailbeats (Fig. 1A), perturbation mode (PM) that consists of high-frequency and low-amplitude oscillations (Fig. 1A), and accumulated mode (AM), which represents the combination of sinusoidal wave (BM) and perturbations (PM) (Fig. 1B): where subscripts bm, pm, and am refer to the parameters in basic mode (BM), perturbation mode (PM), and accumulated mode (AM), respectively.Θ bm and Θ pm are the foil's oscillation angle amplitude in the BM and PM, respectively.f bm and f pm refer to the oscillating frequency in the BM and PM, respectively, and t is the instantaneous time.To control the frequency and amplitude of the perturbation, we introduce dimensionless frequency f = Experiments were performed with a robotic fish swimming in a flow tank that was tethered to a high-resolution force sensor.The robotic fish consisted of a high-fidelity goldfish-like body that was capable of replicating various tailbeat kinematics (see Materials and methods; Fig. 1C).We established the basic tailbeat mode A and B) Foil pitching angle variation in time for the periodic "basic mode" (BM), "perturbation mode" PM (controlled by ( f , Ã) = (10, 0.04)), and "accumulated mode" AM (controlled by ( f , Ã) = (10, 0.04)), respectively, where the time is normalized by 1/f bm .C) Sketch of the experiments in the flow tank with a high-fidelity robotic fish.The motion of the robotic fish's body wave was manipulated with and without perturbations (BM and AM), and the power cost and net thrust were measured.D) A tear-like foil is used to simplify the fish.The foil is of a length L and a head thickness D (diameter of the semicircle), where the pitching center is located at the center of the semicircle.
with a frequency of f bm = 1 Hz and amplitude of Θ bm = 0.349 rad, resulting in a Strouhal number of 0.32, similar to the Strouhal number employed by swimming fish (24)(25)(26)(27).The corresponding Reynolds number Re = UL/ν is 1,900, where L is the body length and ν is the kinematic viscosity of the fluid at a temperature of 20 • .In the experiments, we measured time-averaged net thrust ̅ T, time-averaged power cost ̅ P, and swimming (Froude) efficiency η (= ̅ TU/ ̅ P) of the robotic fish in the BM and AM modes.The power cost, generated thrust, and efficiency ratios with perturbation to those without perturbation are defined as P = ( ̅ P am − ̅ P bm )/ ̅ P bm dimensionless power, to quantitatively compare the performances under different modes.
Experiments were conducted with a perturbation of ( f , Ã) = (4, 0.3).We find P = 265%, T = 366%, and η = 28% (Fig. 2A-C; Movie S1), suggesting that introducing relatively high-frequency and low-amplitude perturbations to tailbeats can improve both net thrust and swimming efficiency simultaneously.To systematically explore how perturbation influences the fish propulsion, we further studied the function of perturbations with f ranging from 4 to 6 and Ã ranging from 0.1 to 0.3.As shown in Fig. 3A and B(i), both P and T are positively correlated with f and/or Ã.This reveals that the stronger the perturbations, the higher the power costs, and the higher the generated net thrusts.Moreover, the normalized net thrust T are all above 0, indicating that perturbations always improve net thrust.In contrast, Fig. 3C(i) indicates that the normalized efficiency η has a complex relationship as a function of f and Ã.In general, with fixed f ( Ã), an improvement of Ã ( f ) will first improve and then inhibit η.

Simulations on pitching foil
To explore the parameter space beyond the range of experimental capabilities, especially with respect to large variations in perturbations ( f and Ã) that would impose high torques on the actuators, we established a computational fluid dynamics simulation environment using FLUENT (ANSYS version 14.0).The simulations also enable a more comprehensive analysis of the flow profiles that result from superimposed perturbations and thus shed light on the hydrodynamic features underlying the observed trends.For simulations, we employed a simplified tear-like foil with the dimensionless thickness of D/L = 0.196 (D denotes foil thickness, Fig. 1D).
To validate the simulation model with respect to physical experiments, we conducted an investigation of the swimming performance of a pitching foil with a perturbation of ( f , Ã) = (10, 0.04) and a basic wave with frequency f bm = 2 Hz and amplitude Θ bm = 0.175 rad at a Reynolds number of 10 3 .In our simulations, we also observed that adding perturbations enhances both thrust and efficiency ( P = 27%, T = 101%, and η = 58%), consistent with our experimental findings (Fig. 2D-F; Movie S2).
To systematically explore the impact of perturbations over a wider range of parameters, we conducted simulations by varying f ranging from 4 to 30, Ã ranging from 0.04 to 0.4, for Reynolds numbers of 1,000, 2,500, and 5,000.Figure 3A-C(ii-iv) summarizes the results after over 80,000 core hours of simulations.We systematically explored the normalized power cost, net thrust, and swimming efficiency as a function of normalized frequency f and amplitude Ã.With increasing in f and/or Ã, P, and T gradually increase, while η first increases and then decreases.As expected, the numerical results conform to those obtained from the experiments, with deviations occurring due to features such as the difference in the profile of baseline tailbeat kinematics, differences in morphology, and differences in dimensions.The similarity in the trends of the improvement in efficiency noticed in both simulations and experiments suggests that the benefits of perturbations are robust and, therefore, could be generic.The optimal case of simulations is obtained at ( f , Ã, Re) = (4, 0.25, 10 3 ), resulting in ( P, T, η) = (196%, 467%, 92%).Additionally, optimal lines for η are formed as Ã ∼ f −1 , indicating that significant improvement in η could be achieved through high-frequency, low-amplitude perturbations and vice versa.Moreover, our numerical results also show that while all perturbations enhance thrust, not all perturbations improve efficiency.Using the power fitting method, we found that the lines formed as Ã ∼ f −1.5 (coefficient of determination R 2 > 0.98) distinguish between improvement and decline in efficiency.

Wake structures and moment analysis
To understand the interaction between pitching foil and flow induced by superimposed high-frequency and low-amplitude perturbations, we further studied the wake structures provided by both BM and AM to visualize the complex dynamics of fluids.At tf bm = 18, we illustrated the contour of vorticity ω * = ωD/U using the numerical simulations results at ( f , Ã, Re) = (10, 0.04, 10 3 ), where ω is the spanwise vorticity.As shown in Fig. 4A, the basic mode (BM) generates 2S wakes (39), a typical reverse Bénard-von Kármán (rBvK) vortex street, leading to a jet-flow (Fig. S2A).In contrast, accumulated mode (AM) generated additional smaller eddies due to superimposing the perturbations (Fig. 4B), resulting in a deflected flow which slants upward first and then downward (Fig. S2B).
There are two stages during the vortex shedding: vortex development (VD) and vortex transportation (VT).VD, viewed as a precursor to VT, physically describes the process of vorticities feeding from the shear layer and is terminated when the generated vortex  2. A-C) Dimensionless time-averaged power, thrust, and efficiency of the robotic fish actuated in the BM and AM mode, respectively, where ( f , Ã, Re) = (4, 0.3, 1,900).D-F) Dimensionless time-averaged power, thrust, and efficiency of the numerical foil actuated in the BM and AM mode, respectively, where ( f , Ã, Re) = (10, 0.04, 1,000).For both experimental and numerical data, the quantities are normalized by the value produced by the BM mode.

Chao et al. | 3
can no longer entrain additional vorticity from the vortex generator (referred to as the vortex pinchoff) (40,41).The vortex pinchoff occurs when the vortex possesses maximum energy from a vortex generator (42), or the vortex generator suddenly stops the energy injection into the vortex, such as the flapping oscillation considered in present work (43).The superimposed perturbations cause the foil to alter its flapping direction rapidly, and multiple vortex pinchoffs occur.As a result, the small-scale eddies provided by perturbations in the wake will be separated between the large-scale vortices generated by the basic oscillation, leading to variations in the positioning of dominant vortices (Fig. S3).
Since the time-averaged power cost can be calculated as dt (see Materials and methods for details), where M f (t) is the torque induced by fluid dynamics, it is worth understanding how M f (t) is affected by superimposed perturbations to demystify the mechanisms underlying efficiency improvements.Figure 5A reveals the M f (t) provided by BM remains the same frequency with θ(t), while the perturbation dominates M f (t) profile in AM.Particularly, a phase lag ϕ exists between θ(t) and M f (t).This ϕ described the phase lag between foil's oscillation and motion-induced fluid dynamics response is actually sourced from the "memory effect" in flow (44).Using the Hilbert transform, we calculated ϕ bm and ϕ am referred to BM-based and AM-based phase lag, respectively, and found |ϕ am | < |ϕ bm |.This discrepancy is attributed to the high-frequency components introduced by superimposed perturbations, which induce rapid transitions in pitching velocity, accelerate the diffusion-convection process in the vorticity layer created on the foil's surface after a change in velocity, and thus facilitate an efficient response of the fluid dynamics to the pitching oscillation.
Mathematically, the minimal ̅ P will be achieved when M f (t) and θ(t) are out-of-phase, which corresponds to M f (t) and θ(t) being in-phase, i.e. ϕ = 0. On the other hand, a smaller phase lag between fluid and flow-driven motion is thought to be beneficial for enhancing the energy harvesting efficiency in a flow-induced vibration system (45).Figure 5B implies that ϕ am is a function of f and Ã and the minimal |ϕ am | values agree well with the data-driven optimal η line in Fig. 3C(ii).This result reveals that phase lag is a critical parameter in affecting efficiency production, a smaller ϕ would cause less energy consumption and potentially improve swimming efficiency.

Analytical models
To explore the mechanisms underlying the improvements in propulsive thrust and swimming efficiency, we further presented analytical models to examine the correlation between dimensionless thrust and efficiency, and the dimensionless frequency and amplitude.
The time-averaged thrust acting on the foil ̅ T is the sum of the time-averaged pressure-induced thrust ̅ T P (> 0) and timeaveraged viscous drag ̅ T V (< 0), as Inspired by previous works on pitching foils (31,46,47), the timeaveraged pressure-induced thrusts read as where S f denotes the surface area of the foil, while can be driven from the Bone-Lighthill boundary-layer thinning hypothesis (48,49).Equation 5 reveals the magnitude of viscous drag produced by the AM is slightly larger than that produced by the BM ( ̅ S4), agreeing with the theoretical prediction.The increase of viscous drag undergoing the AM is owed to the compression of the shear layer and the acceleration of the tangential component of the potential outer flow near the foil.With increasing in Re, ̅ T V approaches zero (Eq.5), suggesting the pressure-induced thrust dominates the thrust generation (50,51).We thus have when the viscous effect is neglected (Re → ∞).As the perturbations introduce the dimensionless amplitude and frequency, the Strouhal number ratio St also increases ( St > 0), and thus timeaveraged thrust is improved ( T > 0, Fig. 3B).Equation 6further reveals the thrust enhancement in AM is mainly sourced from the increase of pressure-induced thrust (Fig. S4).
We next derive the dimensionless efficiency as a function of the dimensionless frequency and amplitude.To do so, we also need to estimate dimensionless power cost, as where M i (t) = J θ(t) is the inertia torque with the mass moment of inertia of the foil J, the M f (t) is calculated as respectively (46).Then, we have the P scaling as (see Materials and methods for details) when A pm ≈ 0 is caused by the low-amplitude pitching.
As a result, we obtain the efficiency η scaling as (see Materials and methods for details) when Re → ∞ and A pm → 0. Equation 9 reveals a smaller ϕ am would lead to a larger η with a special St and Ã.
The analytical models are verified by comparing to both experimental and simulation results in Fig. 6.The scaling equations degenerate as T = 0, P = 0 and η = 0, respectively, when no perturbations have been added ( St = 0, and ϕ bm = ϕ am ).

Discussion and conclusion
Swimming is energetically expensive and animals and robots need efficient kinematic strategies to increase the thrust produced and maximize their locomotion efficiency.This study investigates the impact of introducing kinematic perturbations to the rhythmic tailbeat of fish-like swimming oscillators.By systematically varying the profile of the perturbation in relation to the kinematics, we find that superimposing perturbations consistently increase thrust but not swimming efficiency.Given that only particular perturbations enhance both thrust and efficiency, our study suggests that previous conflicting results may arise from variations in parameter spaces (37,38).
The mechanism we found for improving efficiency complements previous studies.Prior research has suggested that the optimal swimming efficiency is determined by the fluid drag on fins B Fig. 5. A) Normalized foil pitching acceleration θ(t) and flow-induced torque M f (t) in BM and AM, respectively, where ( f, Ã, Re) = (10, 0.04, 1,000).B) The ϕ am contour in f − Ã plane at Re = 1,000, where the solid line refers to the optimal efficiency improvement line in Fig. 3C(ii) and the circle symbol denotes the location of the minimal |ϕ am | at the specific f .and wings, with a lower drag resulting in a higher efficiency (52).However, in our case, we found that superimposing perturbations leads to higher viscous drag but still improves efficiency (Eq.5; Fig. S4).Another previously proposed way to enhance swimming efficiency is to increase the maximal angle of attack (AoA) below the dynamic stall limit, which prevents the shedding of leading-edge vortices (52,53).In our study, perturbations, particularly those with small amplitudes, do not noticeably enhance the maximal AoA (Eq.1); therefore, the efficiency gain from perturbations is not attributable to increased maximal AoA.Additionally, the efficiency improvement may be due to resonance between the natural frequency of the oscillating system and the frequency of vortex shedding (54).We cannot reject the hypothesis corresponding to the resonance effects since reducing the phase lag between body motion and fluid dynamics responses also affects the frequency of vortex shedding.
According to Granger and Paidoussis (44), the phase lag is likely caused by the flow retardation effect, which results in a reorganization of the flow around the body, leading to a time delay between body motion and fluid dynamics response, known as the phase lag.This is different from the time delay caused by the Wagner effect (55) (also known as the "start vortex effect") as the phase lag maintains after the initialization.As the Reynolds number increases, the boundary layer becomes thinner, which in turn diminishes the rate of acceleration for the diffusion-convection of the vorticity.This leads to a smaller reduction of phase lag and thus lower efficiency improvement.Therefore, the higher the Reynolds number, the smaller the improvement in efficiency.Our findings in Fig. 6C are consistent with this, as we observe that the efficiency improvement from certain superimposed perturbations becomes insignificant at higher Reynolds numbers.On the other hand, as Re increases, it would result in a reduction in efficiency improvement since the swimming efficiency provided by the BM increases with increasing in Re (50).This further indicates that superimposed perturbations to improve locomotion may be more beneficial for larval fish than adult fish, since the viscosity has a greater impact on larvae (56).These findings could also explain why larval fish exhibit more random body waves while adults exhibit clearer stereotypical sinusoidal tailbeat kinematics (57).
Further examination of the impact of more general perturbations on thrust and efficiency will be beneficial.Specifically, investigating the effects of nonkinematic higher frequency and lower amplitude random perturbations, which more accurately replicate the perturbations that animals experience in the wild, such as those caused by turbulence during flight or swimming, will be useful.The findings and analytical models presented here lay the foundation for extending to practical applications that can effectively introduce kinematic perturbation to deliver thrust and efficiency improvements to general underwater robots or gliders.

Experimental setup and procedure
Figure 1C shows a schematic of the experimental setup.The experiments were conducted at the Max Planck Institute of Animal Behavior using a flow tank with dimensions of 250 × 875 × 250 millimeters (mm) in width, length, and depth.The flow speed was calibrated to be linearly proportional to the control voltage.Our robotic fish models were designed according to the morphology of the goldfish with a body length of 15 cm.We first scanned real fish to obtain the 3D fish body as a point cloud map.Reverse modeling of the mesh from the cloud map was utilized for our mechanical design using SolidWorks, and the fish body was printed with a 3D printer (Ultimaker S3).Limited by the size of the real fish, we included a single joint in this robot design.The oscillation of this joint is controlled by a waterproof servomotor (Hitec HS-5086WP) driven by a central pattern generator control (58).
Thrust T(t) was measured using a force balance with a fullscale range of 5 kg (Z6FD1, Deutschland HBM).Power P(t) was calculated as the input energy of the motion-generated motor.Data were sampled at a rate of 5,000 points per second using a current acquisition system (NI-9227, National Instruments).Each experimental trial consisted of a 5-s initialization period, during which the robotic fish swam to establish the vortex stream, followed by a 5-s data collection period.To ensure reliability, each set of parameters was repeated five times.For the base mode, we find P bm = 0.775 W, T bm = 0.183 N, η bm = T bm /P bm = 0.236 N/W.

Simulations
In our study, we adopted the mesh system and numerical method previously used by Muhammad et al. (35).The computational domain used in the simulations is shown in Fig. S6.The foil's geometry was defined using the same parameters as in Fig. 1D.The computational domain is a 2D rectangular region with dimensions of 200D × 100D, where D is the foil thickness.The inlet velocity boundary condition is located at a distance of 50D from the pitching center of the foil, with a velocity vector (U, 0).The outflow boundary is situated 150D downstream, with no-stress outflow boundary conditions.The upper and lower boundaries are slip walls located 50D apart from the foil.The computational domain is divided into three grid zones: Zone 1, Zone 2, and Zone 3. Zone 1, inside the blue circle region, has a high grid resolution around the foil.Zone 1 is given an O-xy mesh system where the grids move with the foil as a rigid body.In Zone 1, the first cell is placed at a distance of 0.004D from the foil surface.Zone 2, the black region, has a high resolution to capture large velocity gradients.The blue circle is a circular nonconformal sliding grid interface.It is the interface between Zone 1 and Zone 2. Zone 3, located away from the wake region, has a medium resolution.Zones 2 and 3 employ structured meshes.The prescribed motion of Eq. 1 is controlled using a user defined function with the DEFINE CG MOTION function.
The grid and time-step independence tests were performed to validate the numerical simulations at the most strenuous parameters ( f , Ã) = (30, 0.40) with Reynolds number of Re = 5 × 10 3 .Figure S7 shows the results.The grid independence test is done for three sets of grids with 1.14 × 10 5 (G1), 2.13 × 10 5 (G2), and 4.18 × 10 5 (G3) corresponding to 233, 350, and 525 points on the foil surface, respectively.The time-step of Δt = 1/(2,000f pm ) in the grid independence test is decided by Muhammad et al. (35) and the empirical formula (Δt = min {1/(2,000f pm ), L/‖V‖}, where ‖V‖ is the maximal instantaneous convective flux velocity in the computational domain) proposed by Kinsey and Dumas (59).The dimensionless instantaneous thrust (by ̅ T bm ) over one pitching period (1/f bm ) for the three grid systems show no obvious difference is observed between G2 and G3 (Fig. S7A and Table S1).The difference in the dimensionless time-averaged thrust (by ̅ T bm ) between G2 and G3 is reasonably small, as 0.65%.Mesh G2 is thus adopted.With G2, three time-steps Δt 1 = 1/(1,000f pm ), Δt 2 = 1/(2,000f pm ), and Δt 3 = 1/(4,000f pm ) are tested.The result of dimensionalized instantaneous thrust reveals the difference in the dimensionless time-averaged thrust between Δt 2 and Δt 3 is only 0.53% (Fig. S7B and Table S2).Considering the accuracy and computational resources, Δt 2 is chosen for the extensive work.

Models
In Eq. 3, ̅ T P reads as Here, S f denotes the surface area of the foil, ε(t) is the angular velocity of the fluid around the foil.As introduced before, a phase lag ϕ exits between θ(t) and ϵ(t), we thus have corresponding to BM and AM, respectively.Considering a small Θ bm and even smaller Θ pm , we get cos[θ bm (t)] ≈ cos[θ am (t)] ≈ 1 (Fig. S8), and Θ pm /Θ am ≈ A pm /A bm = Ã.Therefore, the timeaveraged pressure thrusts read as Eq. 4.
The time-averaged viscous drag can be estimated using results of the works on the Bone-Lighthill boundary-layer thinning hypothesis (48,49), as Eq. 5.
We thus have with a positive constant c 0 .Particularly, when the viscous effect is neglected (Re → ∞), Eq. 12 degenerates as Eq. 6.
We next derive the dimensionless efficiency as a function of the dimensionless frequency and amplitude.The time-averaged power cost is considered as the product of the net torque M(t) = M i (t) − M f (t) and θ(t) on the pitching foil, i.e. as introduced before, we get Eq. 7.
For BM and AM, the net torque M(t) can be calculated as Note A pm ≈ 0 caused by the low-amplitude pitching, the Eq. 17 can therefore be simplified as Eq. 8.
As a result, we obtain the efficiency η scaling as In particular, when Re → ∞ and A pm → 0, we get Eq.9.
Chao et al. | 7 see Materials and methods for details), respectively.The Strouhal number ratio is St = St pm /St bm , where St bm = 2f bm A bm /U and St pm = 2f pm A pm /U represent the BM-based and PM-based Strouhal numbers, respectively.
Fig.1.A and B) Foil pitching angle variation in time for the periodic "basic mode" (BM), "perturbation mode" PM (controlled by ( f , Ã) = (10, 0.04)), and "accumulated mode" AM (controlled by ( f , Ã) = (10, 0.04)), respectively, where the time is normalized by 1/f bm .C) Sketch of the experiments in the flow tank with a high-fidelity robotic fish.The motion of the robotic fish's body wave was manipulated with and without perturbations (BM and AM), and the power cost and net thrust were measured.D) A tear-like foil is used to simplify the fish.The foil is of a length L and a head thickness D (diameter of the semicircle), where the pitching center is located at the center of the semicircle.
Fig. 2. A-C) Dimensionless time-averaged power, thrust, and efficiency of the robotic fish actuated in the BM and AM mode, respectively, where ( f , Ã, Re) = (4, 0.3, 1,900).D-F) Dimensionless time-averaged power, thrust, and efficiency of the numerical foil actuated in the BM and AM mode, respectively, where ( f , Ã, Re) = (10, 0.04, 1,000).For both experimental and numerical data, the quantities are normalized by the value produced by the BM mode.

Fig. 6 .
Fig. 6.Scaling for the dimensionless A) power P, B) thrust T, and C) efficiency η.Dashed lines refer to the scaling lines derived from Eqs. 8, 6, and 9, respectively.